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Abstract 

We present an effective theory describing the low-energy properties of an in- 
teracting 2D electron gas at large non- integer filling factors v 3> 1. Assuming 
that the interaction is sufficiently weak, r s < 1, we integrate out all the fast 
degrees of freedom, and derive the effective Hamiltonian acting in the Fock 
space of the partially filled Landau level only. This theory enables us to find 
two energy scales controlling the electron dynamics at energies less than fiuj c . 
The first energy scale, (%uj c /i>) In (vr s ), appears in the one electron spectral 
density as the width of a pseudogap. The second scale, r s hu> c , is paramet- 
rically larger; it characterizes the exchange-enhanced spin splitting and the 
thermodynamic density of states. 

PACS numbers: 73.20.Dx, 73.40.Hm, 73.40.Gk 
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I. INTRODUCTION 



Since the discovery of the fractional quantum Hall effect^ the properties of two- 
dimensional (2D) interacting electron systems in a strong magnetic field have attracted 
persistent attention. Historically, almost all the efforts were directed toward the study of 
the strong magnetic field case when only the lowest Landau level is occupied. Some attempts 
were undertaken! to extend the analysis to a larger filling factor v > 1. However, they were 
all limited to the case when the energy of the electron-electron interaction is much smaller 
than the inter-Landau level spacing, huo c ^> e 2 / (k£), with £ and k being the magnetic length 
and the background dielectric constant respectively. Under this condition, one can describe 
the system in terms of one Landau level only; the effect of mixing of the other Landau levels 
was either neglected or taken into account by perturbation theory!. 

In a weak magnetic field, the typical Coulomb energy exceeds the cyclotron energy huj c . 
In this case, one can start from the Fermi liquid theory in zero magnetic field. The concept 
of Fermi liquid usually enables one to neglect the interactions between quasiparticles when 
determining their energy spectrum in the vicinity of the Fermi level. This concept is based, 
first, on the screening of Coulomb interaction, and, second, on the constraint on the phase 
space allowed for a scattering event!. Landau quantization of the non-interacting quasipar- 
ticles in two dimensions creates a system of discrete levels. The existence of energy gaps 
between the levels affects adversely the screening, which involves the low-energy excitations 
of the electron system. In addition, if the Fermi level coincides with a Landau level, the 
system of non-interacting quasiparticles becomes infinitely degenerate, and even a weak in- 
teraction lifts this degeneracy. This may result in a non-trivial fine structure of the partially 
filled Landau level (PFLL). 

In this parier, we study the low-energy properties of the system with a partially filled high 
Landau levels (the level index N 1 equals to the integer part of [^/2]). By integrating out 
all the other degrees of freedom, we derive the effective interaction between the electrons 
occupying this level. This procedure is valid for a weakly interacting electron gas r s < 1, 
and at sufficiently large filling factors, Nr s ^> 1. (Here r s = \/2e 2 / nhvp is the conventional 
parameter characterizing the interaction strength.) 

The effective Hamiltonian enables us to develop a microscopic theory of the thermody- 
namic and tunnelling densities of states, and a description of spin excitations. In the tun- 
nelling density of states, we find a pseudo-gap with a characteristic width (hco c /2N) ln(iW s ), 
which confirms the result of the hydrodynamic approach!. The thermodynamic density of 
states and spin excitations are characterized by the energy scale r s hu c . Both energies are 
smaller than hu c . It makes our theory complementary to the conventional Fermi liquid 
picture! which is valid in the energy range E Tiu c . 

The paper is arranged as follows. Sec. [H] is a qualitative discussion of the energy scales 
relevant for the problem, and the summary of our main results. Sec. fill] is devoted to 
the rigorous derivation of the effective Hamiltonian describing the low-energy physics. In 
Sees. [TV] and [V] we apply the effective Hamiltonian to study the thermodynamic and the 
tunnelling density of states. The spin excitations are discussed in Sec. jVTf 
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II. QUALITATIVE DISCUSSION AND RESULTS 



Let us first consider an incompressible electron liquid with an integer filling factor v = 
2N. Because of the gap Tluj c in the spectrum, an external in-plane electric field £ can not 
be screened by the 2D electron system. Instead, it causes only a finite polarization per 
unit area, V = \£- The polarizability of the incompressible 2D electron liquid reduces the 
interaction U(r) between two point charges embedded into it, 

U(r)= I ^ ^ —e** (2 1) 

U[r) J (277)2 Kq £{q) e > y ZA ) 

where the two-dimensional dielectric function e{q) is related to the polarizability \ D Y^ 

e(q) = l + ^ x (q). (2.2) 

K 

At small wave vectors, the matrix elements of the dipolar moment between adjacent 
Landau levels, d NtN _ 1: give the main contribution to the polarizability, x ~ n L\dN,N-i\ 2 /^c] 
here ni = l/2~-£ 2 is the electron density on a Landau level. The characteristic spreading of 
the electron wave function on a high Landau level is equal to the cyclotron radius R c = Vp/u c , 
and the estimate of the dipolar moment is cIn,n-i ~ gR c . Substituting \ into Eq. (|2.2j), we 
find: 

e{q) = 1 + $.q, (2.3) 
a B 

where as = h 2 K,/me 2 is the Bohr radius. Eq. (|2.3|) is valid only at small wave vectors, 
qR c < 1. In the opposite limit, qR c ^> 1, a large number of Landau levels participate 
in polarization of the 2D electron liquid. Therefore, the standard Thomas-Fermi screening 
holds: 

e(q) = 1 + —. (2.4) 
a B q 

Formulas ( |2.3| ) and ( |2.4| ) match at qR c ~ 1; the corresponding value of the dielectric constant, 
e(q ~ l/R c ) ~ Rc/ub, is large in the weak magnetic field limit, where R c ^> a B - The 
dielectric function for arbitrary q is given by Eq. ( |3.2(j| ), see Sec. HI . 

As it follows from Eq. (p.3| ), polarization is irrelevant only for interactions on a very large 
length scale, r ^> R 2 /a B , where U(r) is given by the unscreened Coulomb interaction, see 
Fig. [I]. At a smaller scale, R c r R 2 /a B , polarization is important, and Eqs. ( |2.1|) and 
1 yield: 

U(r)^J£). (2.5, 



A t flg < r <C R c , according to Eq. (|2.4j ), Thomas- Fermi screening takes placer, and the 
potential has the form 



2^2 



era 



U(r) ~ —f. (2.6) 
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The renormalized potential Eq. ( |2.1|) is significantly smaller than %uj c . Therefore, unlike 
the bare Coulomb potential, the renormalized interaction does not mix Landau levels. This 
observation enables us to construct an effective theory of the low energy properties of the 
electron system. Similar to the strong magnetic field case, v < 1, the corresponding Hamil- 
tonian is just the energy of interaction between the electrons restricted to a single, partially 
filled Landau level. The main difference is that in our case the interaction potential (|2.1|) 
is much smaller than bare Coulomb potential, and the PFLL has a high index N ^> 1. Al- 
though still unsolvable, this low-energy theory is much simpler than the original one, where 
all the Landau levels were relevant. 

Armed with the effective theory, we are able to estimate the characteristic energy scales 
which control the low-energy dynamics. In order to make these estimates, we employ the 
Hartree-Fock trial function, based on the one-electron wave functions of the iV-th Landau 
level. (See Sees. |IV| and [V| for more details.) At small electron concentration on the PFLL, 
n>e ;$ 1/-Rc' the separation between electrons on this level exceeds the size of the wave 
function, R c . In this case the system under consideration is equivalent to a classical crystal 
of point charges interacting via potential (|2.1|). In the extremely dilute limit, <C a 2 B /R^, 
the polarization is not important, and the usual Coulomb repulsion between the electrons 
of the PFLL is not modified. In this limit the electron system of the PFLL does not differ 
from a Wigner crystal of electrons of the lowest Landau level in the strong magnetic field 
regime. However, due to the polarization effect, there is a wide region of electron densities 
in the weak field regime, a 2 B / ' R\ ^ < 1/-Rc, where the Coulomb repulsion in the crystal is 
replaced by the logarithmic potential ( [2.5|) . In this case, the static properties of the electron 
crystal are equivalent to those of the vortex lattice in a thin superconducting filma. 

We apply the developed picture to study of the energy dependence of the tunnelling 
density of states. We find a suppression of the tunnelling density of states (one-electron 
spectral density) at energies close to the Fermi level Ep, which is similar to the known 
result J^HU for a small filling, factor, v % 1. In the latter case, despite rather different 
approaches, all the authors£2riL3 find a "gap" region of width A g ~ e 2 n l J 2 / n (here n e = n^ =0 
is the total electron density at v < 1). Such a gap can be easily understood: A g equals 
the Coulomb interaction energy of an extra electron with a "frozen" 2D electron system. If 
there were no relaxation of the system at all, the density of states would vanish at energies 
\E — Ep\ ~ A g . The relaxation processes smear out the threshold in the density of states; 
however, the suppression is still strong at these energies, because the characteristic relaxation 
time exceeds h/A g in a strong magnetic field. For a weak magnetic field, there are two 
stages in the evolution of the system after an electron tunnelled in. On the short time 
scale, t ~ 0U7 1 , the polarization of the medium is formed, and the potential induced by the 



additional electron acquires the form ( |2.1| ). At larger time scales, only the electrons of the 
PFLL can be re-distributed. The latter processes are slow, which causes the suppression of 
the tunnelling density of states. The corresponding gap is determined by the interaction of 
the tunneled electron with its nearest neighbors, A g ~ U(r ~ 1/Jn^). For a broad range 
of the filling factors of the PFLL, a 2 B /R 4 c «/if <C l/R 2 c , the width of the gap is: 



It is worth noticing that A g depends on the interaction strength and on the electron density 
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on the PFLL only logarithmically. 

At larger concentrations, ^ -R^T 2 , the single electron wave functions are still orthogo- 
nal to each other, but the densities corresponding to these wave functions have a significant 
spatial overlap. Because of the strong repulsion between electrons at short distances, one 
could naively expect that the interaction on these distances gives the main contribution to 
the energy of the many-electron state. However, for the fully spin polarized electron system 
of the PFLL, the orbital part of the many-electron wave function is antisymmetric, and 
vanishes whenever two electrons have the same coordinates. This suppression of the wave 
function amplitude compensates the large value of the interaction potential at r < ag. At 
larger distances, a# < r < R c , the interaction fl2.6|) decays rapidly, and the corresponding 
contribution to the Hartree-Fock energy turns out to be smaller than the contribution of 
the long-range interaction ( |2.5| ), see Sees. [IV] and |VJ This results in the saturation of the 
logarithmic growth of the gap A g with the electron density: at > R~ 2 the gap in the 
tunnelling density of states is 

A *-^ ln MO- (2-8) 

This estimate was obtained previously in the framework of the hydrodynamic approach!. 

We turn now to the discussion of the spin-flip excitations. The orbital part of the many- 
electron wave function of such an excited state is no longer antisymmetric. Therefore, the 
energy of such an excitation, A s , is determined not only by the bare Zeeman energy but also 
by the extra interaction energy associated with the change in the structure of the orbital wave 
function. Because the orbital part of the wave function is not antisymmetric, the interaction 
on a short range (r < as) now does contribute to the energy of the state. This corresponds 
to the well known exchange enhancement of the g-factoilli Following Ref. |14| we analyze 
the exactly solvable case of a PFLL which has a completely filled spin-polarized sublevel. 
Flipping the spin of a single electron affects the exchange interaction of this electron with all 
the other electrons of the same Landau level. All the lower Landau levels carry equal number 
of electrons of both spin polarizations, and do not contribute to the energy associated with 
the spin flip. Therefore, only l/2iV-th fraction of the total electron density participates 
in the exchange enhancement, and the corresponding contribution to A s is approximately 
N times smaller than the exchange energy per electron at zero magnetic field. For a 2D 
electron gas with Fermi wave vector /cp, the latter energy is of the order e 2 kp/n, and the 
contribution to A s is of the order e 2 kp/ 'kN '. The rigorous result obtained in Sec. [VT| differs 
from this estimate only by a logarithmic factor. For the effective g-factor defined by the 
relation 

9oS = (2.9) 

nuj c 

we obtain 

r /2 3 / 2 \ 

g cS =£ + ^-ln , (2.10) 



if the filling factor v is odd. Here go is determined by the Zeeman splitting in the absence 
of the interaction, and it is usually small (in GaAs the value of g$ is only 0.029). For the 
even filling factors, the value of A s is determined solely by g^. 
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We see that in the weak magnetic field, the energy scale of charge excitations hu c /N is 
parametrically smaller than the energy scale for spin excitations r s %uj c . This is qualitatively 
different from the situation at low filling factor, v % 1, where both these excitations are 
characterized by the same energy scale e 2 / k£. 

III. EFFECTIVE HAMILTON! AN 

We start from the full Hamiltonian of the system 

H = H + Hi nt , (3.1) 
where the Hamiltonian for non- interacting electrons, Hq, is given by 

H o = J2 i nhuJ c ~ Vi.kV'n.fc- ( 3 - 2 ) 

n,k 

Here n and k are the Landau level index and the guiding center coordinate respectively (we 
adopt the Landau gauge). The fermionic field operators in Eq. (|3.2|) satisfy the standard 
anticommutation relations {ipn.ki^m.^} = 0, V'm,?} = 5 mn 5k q , and p is the chemical 

potential (we include the energy hu c /2 into p). To simplify the notation, we omit the spin 
indices. We will also neglect the Zeeman term goahuo c /2 in all the intermediate calculations. 
This is legitimate because go <C 1 and therefore all the effects associated with the small 
Zeeman energy can be included in the effective Hamiltonian after the fast degrees of freedom 
are integrated out. 

Electron-electron interaction is described by 

Hint = — I I i r : (p(r-i) - n e ) (p(r 2 ) - n e ) : , (3.3) 

2k J J \r 1 — r 2 \ 

where : . . . : stands for the normal ordering, p is the electron density operator, p(r) = 
W(r)ty(r), and n e is the average electron density. The electron annihilation operator ^/ 
is related to the operators ip n ,k by ^(r) = J2 n ,ki J n,k l Pn,k('r), where <p nt k(r) are the single 
electron wave functions in a magnetic field. 

In order to develop the effective low-energy theory, we notice that the expression for the 
partition function, Z = Tr exp(—/3H), can be written as 

Z = Tr N (Tr' N exp(-[3H)), (3.4) 

where Tr^ is the trace over the Fock space of the PFLL (N is the index of the PFLL, f3 
is the inverse temperature) and Tr' N means trace over all the other Landau levels. At low 
temperatures, f3Tiu c 3> 1, Eq. (|3.4j ) acquires the form 

Z = Tr N exp(-0H ef f), (3.5) 

where the effective Hamiltonian H e ff is defined as 

% = -limilnA(A), A(p ) = Tr' N {e-^\. (3.6) 

On— >oo In *■ ' 
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Let us emphasize that, in general, the low frequency properties can be described only 
by means of the effective action which includes a retarded interaction. The description by 
an effective Hamiltonian with instantaneous interaction is possible only if the characteristic 
time of the retardation is much smaller than h/E c , where E c is the maximal energy scale 
which is considered within the effective theory. This condition is met for the problem under 
consideration. Indeed, the summation in Eq. ( |3.6|) involves inter-Landau level transitions 
only. These transitions are associated with the energy gaps not smaller than hu c . Corre- 
spondingly, the characteristic time of the retardation related to such a transition is t ^ ^T 1 . 
On the other hand, we saw at the Sec. [TJ that the largest characteristic energy E c = A s is 
still parametrically smaller than hu c , which enables us to neglect the retardation. 

We now turn to the calculation of the effective Hamiltonian (|3.6| ). The simplest way 
to proceed is to use the Hubbard- Stratonovich transformation: first, we introduce auxiliary 
scalar field <f>(r, r) to decouple quartic in \l/ interaction part of the Hamiltonian; second, we 
perform the summation over all the fermionic degrees of freedom not belonging to PFLL. 
After that, we integrate out auxiliary fields by means of the saddle point approximation and 
thus obtain the effective Hamiltonian H e ff. 

Below we implement this procedure. Performing the first step of the Hubbard- 
Stratonovich transformation, we find for A from Eq. Q3.6| ) 

Tr' N { e -A*»T T exp (- J d 3 ^) (A*)-".)) } » ( 3 - 7 ) 

where T T is the imaginary time ordering, and 3D vector £ represents time and space coor- 
dinates, £ = (t, r). The domain of integration over £ is: r £ [0, /?o], x £ [0, L x ], y £ [0, Ly], 
where L x x L y is the geometrical size of the 2D system. We will omit the Planck constant 
in all the intermediate calculations. 

The measure of the functional integral in Eq. ( |3.7| ) is determined by the condition 

fV(j> exp^l^y'd 3 ^!)^)^!-^))^!, (3-8) 
and the function K is defined by the equation 

/ ^ 3 ^(6-6) e2 f (T3 ~ T2 | ) = 5( (3.9) 
J K\r 3 —r 2 \ 

The density operator is 

p(r,r) = V(r,T + e)V(r,T), e -> +0, (3.10) 

where the fermionic operators in Matsubara representation are ^(r) = e THo ^fe~ T , and 
^(t) = ^ r '(— r). By introducing the infinitesimal positive time shift e into definition ( |3.10| ) 
we preserve the normal ordering of the fermionic operators in the interaction Hamiltonian, 
see Eq. Q). 

After the decoupling is performed, the fermionic part of the Hamiltonian in Eq. ( |3.7| ) 
becomes quadratic, which enables us to make the second step of the Hubbard- Stratonovich 
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transformation, i.e. to carry out the trace Tr' N over the "fast" fermionic degrees of freedom. 
The result depends solely on the variables belonging to the PFLL: 



A(A)) = e~^jV(j> e- Fo ^T T \ e -^M e ~h{^N^ N } 

N 

7T 



(3.11a) 



cf 6 K60(6)0(6) (Kfa-b) -4(^(6,6)^(6,6: 



(3.11b) 

A {0, M = y d 3 6>(£) (p N (Z) - <) > (3-llc) 
F 2 {0, f ^} = J d 3 ^ J d 3 & [f N (t+)* N fa) - 2G N {£ 2 , 6)] 0(6)^(6, 6)0(6), (3.11d) 
Gjv(6,6) = -PN{r^r 2 )e^-^ N ^e{T X - r 2 - e), (3.11e) 



where the 3D vector £ + is (r + e, r). In Eqs. (|3.11| ) \&jv is the projection of the fermionic 
operator (in Matsubara representation) on the PFLL, ^n{$) = e T ^~ NuJc ' J2k i } N,k i fN,k(i'), 
operator ^jv(r) is defined as ^^(r) = & n (—t). Projections of the density operator p and 
of the Hamiltonian H on the PFLL are given by pn(£) = ^?n(£ + )^n(£) and by Hq = 
(Nu c — p) J d?r_fiK[{r) respectively. In Eq. (|3.11ej ), -P/v(r2,ri) is the projection operator 

n e — N/tc£ 2 is the average electron concentration on the PFLL. 



n 



N 



onto the PFLlJl 

Operator F 2 in Eq. ( |3.11d|) is defined in such a way T T F 2 = F 2 . This follows from the 
anticommutation relation for the fermionic operators j\I/jv(T*i), ^jv( r 2)} = Pn{tx, t 2 ). 
The Green's function G^(6,6) is the kernel of the integral operator: 



Ga, — P, 



N 



--Q- +P + ~£ ~ Pn^4>Pn 



Ni 



(3.12) 



Here Ph is the projection operator onto the space of functions orthogonal to the states of the 
PFLLEil, a nd Ti^ is the Hamiltonian of an electron in the magnetic field and in the external 
potential (f>: 



2m 



r,r . 



(3.13) 



The function G^ describes the evolution in the external field <ft of electron states constrained 
to the Landau levels n ^ N . This constraint is implemented by the introduction of the 
projection operators into Eq. (|3.12|) . At n = r 2 the function G^ is defined as 



G^in, n) -> g^ti, ri + e) 



(3.14) 



in accordance with Eq. (|3.10|) . 

The first term in Eq. (|3.11b| ) results from the summation over the fermionic states with 
n 7^ N, and represents the "thermodynamic potential" of electrons in these states in the 
field 0. The function /^(^i, 6) is defined as the kernel of the operator 



2 In 



(3.15) 



S 



with from Eq. (|3.12 ). Here, the factor 2 comes from the spin degeneracy. 

In order to perform the integration over the auxiliary field in Eq. ( |3.11a| ), we employ 
the saddle point approximation. It means that one has to expand the functional i*o{0}, see 
Eq. ( |3.11b|) , and the operator F 2 {0, ^n, ^v}, see Eq. ( p.lldp , up to the second order in 0. 

We start from the expansion of the operator F 2 {0, ^tv, ^tv}. It contains explicitly a 

factor, which is bilinear in 0, and therefore we can replace in Eq. ( |3.11d|) by the Green's 
function Go of an electron in the absence of an external field, = 0. Then, Eq. ( |3.12 ) 
enables us to find the explicit form of Gq\ 



Go(ri,7i;r 2 ,<r 2 ) = £ ^^(v^-^^ 

k,n^N 

x [0(t 2 - n)e(N - n) - 6{n - r 2 )9(n - N)}. 



(3.16) 



When expanding the first term in Fo{0}, we use definition (|3.15|) and the solution of 
Eqs. ( |3.12j ), ( |3.13| ) up to the second order in 0, and we obtain: 



//V r 

+ / d 3 ^ [ d 3 60(£i)0(6)G o (6, £i)G (£i, 6) 



(3-17) 

where Qq = [(N — l)hu c /2 — /i] NL x L y /ir£ 2 is the thermodynamic potential of the system 
of noninteracting electrons on the completely filled Landau levels. The third term in the 
functional i 7 o{0}, see Eq. ( 3.11b| ), is already bilinear in 0. 
Substitution of Eq. ( |3TTTD into Eq. ( pTb| ) yields 



f o {0} = p n Q - i J J d^d^i^s - 6) 0(6), 



(3.18) 



where the second term describes fluctuations of the field which are renormalized due to 
the integrated out degrees of freedom. The kernel S equals to 



Sir, 



duid 2 q nq 
(2tt) 3 2vre 2 



e(q,u)e l 



qr— iu)t 



(3.19) 



where e(q, uS) coincides with the dielectric function calculated in the random phase 
approximation!!! : 



e(q,u) 



1 



Kq 



U(q,u) 



(3.20) 



Here, the polarization operator Tl(q, uj) in Eq. ( p.20|) is given by 
U{q,uj) = 2 f d 2 rdre~ icir+iu}T \G {r, r; 0, 0) + 



G N (r, r; 0, 0) G (0, 0; r, r) + G N (0, 0; r, r] 



(3.21) 



where the functions Go and Gat are defined by Eqs. (|3.16|) and ( |3.11e|) respectively. The 
explicit form and the asymptotic behavior of U(q,uj) are presented in Appendix [A[ 
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Within the approximations made, we can transform Eq. ( |3.11a|) to the form of a Gaussian 
integral over the auxiliary field 0. The result is 



-Fi{<j>,p N }-F°{<j>,^ N ,-<Si N } 



A(f3 Q )=e- Mno+H ^JV<pT T 
x exp Q J d 3 ^ J d 3 60(6)0(6)«5(6 -6)) , 



(3.22) 



where F 2 ° * s obtained from operator F 2 , see Eq. ( p.llclp , by replacing with Go- 

Before we proceed further, let us emphasize that the saddle point approximation can be 
justified only if the characteristic value of the fluctuations of the auxiliary field is small 
enough. The estimate which will be presented later in this section shows that the expansion 
quadratic in is parametrically valid for the weakly interacting electron system at large 
filling factors. 

Now, we are in the position to perform the actual integration in Eq. (|3.22|) . Operator iq, 
see Eq. ( j3.11c|) , is a linear functional of 0, whereas operator F° , see Eq. (|3.11d|) , is quadratic 
in 0; because the fluctuations of field are small, the typical value of F 2 is much smaller 
than that of iq. If we neglect the term F® at all, the functional integration in Eq. ( |3.22| ) 
can be easily performed. The calculation is further simplified if we approximate S(t) ~ 
5(t) J S(r)d,T. It is valid approximation for the description of the low-energy dynamics of 
electrons belonging to PFLL: S is rapidly decaying function of time: S(t) cx exp(— u; c |t|) at 
" 1 , as it follows from Eqs. Q3.19| ) - ( |3.21| ). The integration over field then results in 



> 



OJ. 



an exponential exp (—/3oH^/\ with the Hamiltonian of density- density interaction, 



H 



eff 
hit 



d 2 r 1 d 2 r 2 U(r 1 - r 2 ) 



(3.23) 



x : (piv(rq) - nf) (p N (r 2 ) - nf ) 



The operator ( |3.23|) acts within the Fock space of PFLL, and the renormalized pair inter- 
action U{r) is related to <S(£) by 



/ d 3 ^S(^-^)U{\r 3 -r 2 \)=5{ 7q-r 2 ) 



(3.24) 



Substitution of Eq. ($J% into Eq. ( gg§ ) yields: 

d 2 q 2vre 2 



U(r) 



(2ir) 2 nq e(q) 



3 «qr 



(3.25) 



Here the static dielectric function e(q) = e(q, to = 0) describes the renormalization of the 
bare Coulomb potential due to the integrated out degrees of freedom. This function can be 
easily found with the help of Eqs. (|3.20|) and ( |A7a| ), and it has the form 



e(q) = 1 + 1 - Jl (qR c ) 



qa B 



(3.26) 



for the large filling factor N ^> 1 and for values of wave vector q much smaller than Fermi 
wave vector, hp. In Eq. (|3.26|) , as = h n/me 2 is the Bohr radius and J7q(x) is the zeroth 
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order Bessel function. Asymptotic behavior of the renormalized potential U (r) was discussed 
in Sec. H, see also Fig. [l]. 

The term F® from Eq. ( |3.22| ) can be taken into consideration by means of the perturbation 



theory using the small parameters r s and 1/JV. We show in Appendix |B| that this operator 
generates term of the effective Hamiltonian which is linear in density. This linear term 
can be trivially included into the free electron Hamiltonian H% by a shift of the chemical 
potential. 

Finally, substituting the resulting expression for A in Eq. ( |3.6| ) and restoring the Zeeman 
energy term, we obtain the expression for the effective Hamiltonian: 



H eff = n-^j d 2 rp N (r) + H e H+ (3.27) 
1 r 

-g huj c J d 2 r (*iv a (r),^;(r) - ^ T (r),^ T (r) 



where H^/ is defined by Eq. ( |3.23|) , the explicit expression for the shifted chemical potential 
/i* is presented in Appendix 10, see Eq. (|B30|), and the thermodynamic potential Q of the 



electrons on the filled Landau levels is given by Eq. ( |B26|) . We will not need concrete values 
of fi* and Q in the further calculations. 

Hamiltonian ( |3.27| ) is the main result of this section. The physical meaning of this 



Hamiltonian is that the low-frequency dynamics of the system is described by the electrons 
belonging to the upper Landau level. Interaction between these electrons is renormalized 
due to the large polarizability of all the other Landau levels. 

Let us now discuss the condition of the validity of the approximations we made. As we 
mentioned earlier, the saddle point approximation can be justified only if the characteristic 
value of the fluctuations of the auxiliary field is small enough. Below we show that the 
fluctuations at all the spatial and time scales are parametrically small for a weakly interacting 
electron system at large filling factors. 

The magnitude of the fluctuations localized within the spatial range r and within the 
time interval r can be estimated from Eqs. (|3.22|) and ( P-19|) as 



h ( » 2 



e 



T,r t \ nreir 1 , r v 



(3.28) 



Further estimates depend on the relation between the scale r and the cyclotron radius R c . 

If this spatial scale is large, r 3> R c , only the dipole transitions between the nearest Lan- 
dau levels are induced. The corresponding off-diagonal matrix elements of the Hamiltonian 
are of the order of (R c /r)\4>\. The mixing of the Landau levels by the field <fr is small if 



(hu> c + h/r) 2 



)t ' t ( —\ <1. (3.29) 



Combining Eq. (gjg ) with Eq. (|3~29|) and using the results (|3.20|) and ( |A7cj ) for the dielectric 



function, we find the condition of small fluctuations 

r s N > 1. (3.30a) 
Therefore, the long-range fluctuations are not "dangerous" in the weak magnetic field regime. 
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Let us now analyze the short range fluctuations, r <C R c - In this case transitions between 
distant Landau levels are possible, and the requirement of the smallness of fluctuations 
coincides with the standard one for the validity of the random phase approximation at zero 
magnetic field, 

r s < 1. (3.30b) 

Thus, the saddle point approximation is valid for the weakly interacting electrons in a weak 
magnetic field. 



In the following sections, we will apply the effective Hamiltonian ( 3.27 ) to describe various 
physical effects associated with the PFLL. 



IV. GROUND STATE ENERGY AND THERMODYNAMIC DENSITY OF 

STATES 

In this section, we evaluate the ground state energy and thermodynamic density of states 
dn e /dp as a function of the filling factor of the partially filled Landau level Au = u — 2N. 
We consider explicitly the case Av < 1. System at the filling factors 1 < Au < 2 can be 
analyzed with the help of the electron-hole symmetry, and the the thermodynamic density 
of states for this case can be obtained from the results for Au < 1 by the replacement 
Au ->• 2 - Au. 

We assume that the ground state is spin polarized at Au < 1, and thus omit the irrelevant 
Zeeman term0. The chemical potential /i* for the given electron concentration on the PFLL 

= Mj L x L y is found by differentiating the ground state energy of Hamiltonian (|3.27|) with 
respect to the number of electrons M: 



dE (M) 



01^/10 



= — w-/ E (M) = x - ' /M , (4.1) 
P dM ' ov ; (0|0) M V ; 

where |0) M is the wave function of the ground state of the system with M electrons on the 
PFLL. 

Evidently, it is sufficient to consider only the filling factors Au < 1/2. At larger filling 
factors, 1/2 < Au < 1, one can use the electron-hole transformation within the spin sublevel, 
= \lzjv, and study the system of holes with the filling factor 1 — Au described by the 
Hamiltonian 

H ef f = {l-H*Jd 2 r ( JL - fa(r?) + H*H, (4.2) 



1 
2 



d T\d r 2 U(r 1 -r 2 ) 



p h N (r,)-n N h 



p h N (r 2 )-n N h 



Here p h N = ^% is the hole density operator p 1 ^ and = 1/(2tt£ 2 ) —n^ is the average 

density of holes at the PFLL. The first term in H^t corresponds to the shift of th^ chemical 
potential of the completely filled Landau level due to the exchange interaction^ 
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2tt£ 2 / d 2 rU(r)P N (r,0)P N (0,r). 



(4.3) 

Here Pn^i,^) is the projector operator on the PFLlM The integration in Eq. ( f4.3| ) can 
be performed with the help of the explicit form ( [B21|) for -P/v(r!,r 2 ). Under the conditions 
( |3.30| ), the calculation yields 



A ex = hu c JL In 



(4.4) 



It follows from Eq. ( |4.2| ) that filling of one spin sublevel results in the shift — A ex of the 
chemical potential. This enables us to find the average thermodynamic density of states 
which is defined by the relation 



. dfj. 



1/2tt£ 2 



li(Av = 1-5) - ii(Av = 5) 



+0. 



(4.5) 



With the help of Eq. (O), we obtain 



dn e 
4 dn 



m 



h 



\/2r s In 



' 2 3/2^ 



(4.6) 



It is worth noticing, that the thermodynamic density of states ( |4.6| ) does not depend on 
magnetic field at all. We will see below that for a broad range of filling factors of the PFLL 
Au, the actual value of the thermodynamic density of states is close to the average one given 
by Eq. (Ep|). 

The true ground state of the Hamiltonian ( pj.27| ) for non-integer filling factors is not 
known. In order to estimate the ground state energy, we use the Hartree-Fock trial function 
analogous to that used by Maki and ZotoJll for the lowest Landau level 



M 
3=1 



■0. 



R 



|0)m 

d 2 r$ K (r)¥ N (r) 



o ■ 



(4.7) 
(4.8) 



where $r(t*) is the normalized one-electron wave function of a coherent stateE! on the PFLL 
with the guiding center localized about point R. In the Landau gauge, A y = —Bx, A x = 0, 
this function has the form 

3(r) = ^{^) N e^\ (4.9) 
rapidly decreases with the distance between 



The overlap between two functions 
their guiding centers: 



|rf 2 r$^(r)$ R2 (r) = e- 



■i[RixR 2 ]z/2£ 2 e -(R 1 -R 2 ) 2 /« 2 



(4.10) 
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Equation ( [4.10| ) shows that the overlap is exponentially small even when the distance between 
the guiding centers is smaller than R c and the electrons are not separated in space. 

For small filling factors Av < 1/2, we can choose the guiding centers Ri separated by 
the distance much larger than the magnetic length I which enables us to neglect the non- 
orthogonality of the coherent states. Then, the expression for the energy Eq(M) given by 
Eq. (|4.1|) acquires a simple form 



M 



E (M) -- 
Vhf(R) 



-Y^V HF {\R t -R 3 \) 
V H (R)-V F {R), 



Mn 



N 



d 2 rU{r), 



(4.11) 



where Hartree Vh and Fock V F potentials are given by 



V H {R)=JJd 2 rd 2 r' U{r-r') |$ (r)| 2 |$ R (r 
d 2 q 2vre 2 1 



/M2 



(27t) 2 Kq e(q) 



<? 2 £ 2 +iqR 



(4.12) 



and 



V F (R)=JJd 2 rd 2 r / U(r-r')%(r)§* R (r')%(r')<S> Il { 
d 2 q 2vre 2 1 



(27r) 2 nq e(q) 



.2 p2' 



(4.13) 



respectively. Here L^{x) is the Laguerre polynomial^. We used Eqs. ( pj.25| ) and ( |4.9| ), when 
deriving Eqs. ( [4.1 2|) and ( |4.13|) . Small corrections appearing due to the non-orthogonality 
of the functions (fl9|) were discussed in Ref. |18. 



Now, the energy given by Eq. ( [4.11| ) should be minimized with respect to the positions of 
the guiding centers FL. The best configuration corresponds to the guiding centers arranged 
in a triangular lattices. Below we present analytical results for the energy of the variational 



-1/2 



ground state ([4.7|) in different domains of filling factor Av. 

Dilute system, Au ^ N^ 1 - In this case, the lattice constant a = (^y/3n^ /2^j is 
large, a 3> R c , so that the Fock potential ( [4.13|) is exponentially small and the Hartree term 
Vh{R) coincides with the interaction potential U(R). In this limit, two situations may be 
distinguished. 

For an extremely dilute system, Au <^ N~ 3 rj 2 , the lattice constant is larger than R 2 /cib, 
and therefore the Coulomb potential is not renormalized. In this case, the ground state 
energy coincides with that of the Wigner crystal on a neutralizing background, 



E Q (M) = -a 



e 2 M 3 / 2 



(4.14) 



where a = 0.782 ... is a numerical constant^, and M$ = L x L y /2ir£ 2 is the number of states 
on a Landau level. The use of Eq. (|4.1| ) immediately yields for the thermodynamic density 
of states in this regime 



14 



dn e 



0.542 m 



N 



(4.15) 



dji r s h 2 

The negative thermodynamic density of states for an analogous system was considered earlier 



in Ref. 22 



For a moderately dilute system N~ 3 r~ 2 <C Av <C AT -1 , the typical distance between 
electrons on the PFLL is much smaller than R 2 /a B . At such distances the potential of 
interaction is strongly renormalized by the screening and it is given by Eq. fl2.5|) . With the 
logarithmic accuracy the energy of the system is given by 



E (M)~-M^ In 
ov ; 8N 



3„2 



N r 



M 



(4.16) 



and, correspondingly, the thermodynamic density of states for moderately dilute system 
takes the form: 



dn e 
<9/i 



4 TYh . . _ —•. 



(4.17) 



"Dense" limit, N^ 1 <C Av <C 1/2 - In this case the distance between the nearest guiding 
centers is smaller than R c . The asymptotic behavior of the Hartree and Fock potentials 
depends on the relation between the magnetic length £ and the Bohr radius a B . We restrict 
ourselves to the case £ ^> a B , or 



Nr 2 s > 1, 



(4.18) 



which enables us to simplify Eqs. ( [4.12j ) and ( |4.13| ) 

3e 2 a B 



V H (R) 
V f (R) 



e 2 a B 



+ -^ln(^ 



2irR c R 4ir 2 R 2 c 
e 2 a B 



2irR c R(l + Ra B /2£ 2 ) 



R 
R 



+ e2(lB l n (El 
R 2 \a B 



(4.19) 



for the range of distances £ <C R <C R c - 

The resulting potential Vhf{R) = Vh(R) — Vp(R) is a smooth function of R at 
R <C £ 2 /a B . This enables us to find the analytical expression for the ground state en- 



ergy Q4.11p in the region N 1 r s 2 <C Av <C 1/2. First, we can approximate J2j^o VHF(Rj) 



rig / d RVhf(R) — Vhf{o)i as the potential varies only logarithmically over the lattice cell. 
Second, we notice that the spatial average of the Hartree potential is exactly equal to the 
average of the bare potential U(r), and these two terms cancel each other in Eq. ([4.11 ) for 
E (M). Finally, the integral of the Fock potential is proportional to the exchange shift, A ex , 
see Eq. ([4.4|). The resulting energy of the system is 



E (M) 



M 2 



2Mz 



-A, 



(4.20) 



M- 



W 



JL hi ( —n) + ln (Nr s ) + —L- 



which yields the thermodynamic density of states of the form 



15 



dn e 



m 

fc2 



y/2r B In 



' 2 3/2^ 



X 



3V2 



327rr, In 



(4.21) 



The leading term in Eq. Q4.21|) coincides with Eq. ( }4.6| ). 

The filling factors Ai/ > 1/2 can be considered similarly by using the electron-hole 
symmetry, see Eq. (|4.2j ), and the trial function of type (|4.7| ). The overall dependence of the 
thermodynamic density of states on the filling factor is shown in Fig. 0. 



V. SPECTRAL DENSITY 



The value of the one electron spectral density A(e) can be measured in the tunnelling 
experiment JHH. For the electron states on the PFLL, A(e) is defined by relations: 



A(e)=A+(e)+A-(e) 



(5.1) 



A' 



M, 
1 

M 



An 



6(E fi 



E, 



- ej 



EK^I^|0)| 2 5(E -^-e) 



* k,; 



Here, |m) is an eigenstate of the full Hamiltonian ( |3.1| ), E^ is the corresponding eigenvalue 
and rh = stands for the ground state of the system. Functions A + (e) and A~(e) describe 
the introduction of an extra electron or an extra hole onto the PFLL respectively. 

We are interested in the spectral density at energies |e| -C hu c . All the behavior of the 
system at such low energies can be described by the effective Hamiltonian ( p. 27 ). Following 



the method of Sec. |TJ, we can express the spectral density (|5.1|) in terms of the eigenstates 
of the effective Hamiltonian as 



A+ ^ = M:J2\{ m An 0) 2 S(E m -E -e 



* k,r 



A-(e) 



M, 



J2\(m\^ k , N \0)\ 2 S(E -E m -e) 



* k, 



(5.2) 



where m denotes the eigenstates of the effective Hamiltonian (|3.27|) and E m is the corre- 
sponding eigenvalue, m = is the ground state. The Z-factor in Eq. (|5.2|) describes the 
overlap between the low-energy eigenstates of the full Hamiltonian and the eigenstates of 
the effective Hamiltonian. The value of Z-factor coincides with the quasiparticle weight for 
a 2D degenerate plasma in the absence of the magnetic field, 

3r, 



Now we use Eq. ( |5.2| ) to evaluate the spectral density. Our analysis is based on the trial 
function (|4.7| ). Even though we are not able to find the details of the energy dependence of 
the spectral density, this approach provides a reliable estimate for the energy scales involved. 
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Because the trial function ( |4.7|) is constructed from the coherent states (|4.9| ), it is con- 
venient to rewrite Eq. ( |5.2|) in terms of the operators iftn, V4t> see Eq. (|4.8| ), creating or 
annihilating an electron in the coherent state: 



ZTT 

g + (R,t) = (o 



S- I dte- id 



-itH, 



d 2 R 

L X Ly 



eff tpn.e eff ii 



R 



JtH, 



(5.3a) 

(5.3b) 
(5.3c) 



Here we used the representation of the projection operator in terms of the coherent states 



PN(r 1 ,r 2 ) 



d 2 R 



$R(r 2 )$R(ri). 



Let us concentrate below on the calculation of the spectral density for the hole excitations 
A~(e). (The calculation procedure for A + (e) is similar and will be briefly outlined later.) 
At small time t, expression (|5.3cj) can be approximated by the formula 



g-(R,t) « 0-(R,O)exp(z*£T(R) - ]- [tD(R)] 2 



(5.4) 



with parameters 



E~(R) 



[D(R)f 





A, 








(0 


^rV'R 





(5.5) 







H eff, 


Heff,1pR_ 






(0 


V'r.^R 





Substitution of Eq. 



into Eq. ( |5.3a| ) gives the result 



Z fd 2 R g-(R,o) 



2tt J L x L y D(R) 



-exp 



(S-(R) 



2 [D(R)]< 



(5.6) 



for the spectral density. 

Now we use Maki-Zotos trial state ([4.71) to find the functions in the integrand in Eq. ( |5.6D . 
At a small filling of the PFLL, Au <C 1, we obtain with the help of Eq. ( |4.10| ) 



^-(R,0)=^e-( R - R ^ 2 ^ 2 
r, 



(5.7) 



which describes the "density of the centers of orbits" in the Wigner crystal. Energy E~(R) 
defined in Eq. ( |5.5| ) turns out to be independent on the position R because all the electrons 
in the Wigner crystal have the same energy: 



E~ = £ V HF (\Ri 

R^O 



nf J d 2 rU{r) - fi*. 



(5.8) 
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Here the Hartree - Fock energy Vhf is denned by Eq. (|4.11|) , and the second term describes 
the neutralizing background. The parameter -D(R) is also independent of R and it is given 
by 

D 2 = £ 2 J2 [V^F(|Ri|)] 2 . (5.9) 



Because none of the parameters in the exponent in Eq. (|5.6|) depend on R, the integration 
can be easily performed, and it yields: 

^ A ^ ex +^}' (510) 

Let us now discuss the physical meaning of the energies E~ , D. After an electron is 
removed from the PFLL, the system acquires extra energy equal to — E~(R.) due to the 
interaction between the hole and all the electrons of the PFLL. If all the other electrons of 
the PFLL were "frozen", this state would give a contribution oc 5(e — E~ (R,)) to the spectral 
density. (Classical model of the frozen electrons was used in Ref. [II] for the electrons on 
the lowest Landau level.) However, the state formed right after the tunnelling is not an 
eigenfunction. The decay of this initial state leads to the finite quantum width D of the 
spectral density, see Fig. |^ a. 

To find the region of applicability of expression ( |5.10|) we notice that the expansion in t we 
used in Eq. ( |5.4j) is valid for \t\ <C \E~ / D 2 \. The characteristic time determining the value of 



the integral in Eq. (5.3a) can be estimated as t ~ |e — E \/D . Combining these two condi 



tions, we find that Eq. ( |5.10|) is applicable for \e — E | <C \E |. At the boundary of applica- 



bility, the spectral density ( |5.10| ) is proportional to exp [— (E /D) 2 ], and therefore Eq. ( |5.10| ) 



describes the main contribution to the spectral density only if \E~\ ^> D. Substituting the 
explicit expression for Hartree- Fock potential ( |4.11| ) into Eqs. (|5.8|) and (|5.9|), we see that 
the ratio D/\E~\ is maximal at Au ~ 1/2, where its value is D/\E~\ ~ 0.1(ln N)^ 1 / 2 <C 1. 
Furthermore, we believe that the use of the Hartree-Fock function sets the upper limit for D. 
This is because trial function ([4.7| ) does not take into account the correlations in the motion 
of centers of orbits, and thus overestimates zero-point fluctuations of the electron density. 
The use of a more sophisticated!! trial function, or of a_phenomenological model that has 
the correct spectrum of excitations of the Wigner crystall!3 would restrict the summation in 
Eq. ( |5.9| ) to the nearest neighbors only, and further reduce the value of D. 

The spectral density of holes has a sharp peak, because a hole corresponds to a vacancy 
in the Wigner crystal. In contrast, the tunnelling electron may create an interstitial at an 
arbitrary position in the elementary cell. This gives rise to a finite width of the spectral 
density A + (e), even for the "frozen", D = 0, crystal. The energy band for an interstitial is 
given by 



£ + (R) = E 



R; - 



V HF (\K l -R\) 

l_ e -(n-R l ) 2 /2fi 



nf \d 2 rU{r)-^. (5.11) 



The leading term in the quantum width D for electrons does not depend on R, and it is 
still given by Eq. (p^9|). 

To calculate the spectral density A + (e), one can use Eq. (|5.6| ) with E~(R) and {?~(R, 0) 
replaced by E + (R) and 1 — £/~(R, 0) respectively. The width of the electron spectral density 
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is determined by the energy band ( |5.11|) . The abrupt edge of the spectral density at e = 
minR {E + (R,)} is smeared only by a small quantum width D, see Fig. ||] a. 

The minimal energy necessary to create an interstitial, min{£' + (R)}, is of the order of 
energy \E~\ needed to create a vacancy. Therefore the difference A g = min{.E + (R)} — E~ 
exceeds significantly D and thus determines the pseudogap in the spectral density fl5.1|). 
Using Eqs. (|5.8|) , ( |5.11|) , and the asymptotic behavior of the Hartree-Fock potential Vhf(R), 
see Sec. [TV], we obtain the width of the pseudogap for different filling factors: 



A, 



IN 

lUJc 

< 2N 



In 



N 3 r 2 Av 



In (Nr s ), 



Av < N~ 3 r~ 2 
jy~3 r -2 ^ Az/ < N~ l 

N" 1 < Au < 1/2, 



(5.12) 



The leading term in the width of the gap agrees with the result of Ref. |6| for the broad range 
of filling factors Au ^ 1/N. The overall dependence of the characteristic energy scales on 
the filling factor Av is shown in Fig. |3] b. 

It is important to notice that in the dominant region of filling factors, Av > 1/N, the 
width of the gap A g is determined by the renormalized interaction ( |2.5| ) at distances r ~ R c , 
which are large compared to the characteristic distance between the electrons of PFLL. The 
gap is insensitive to the short range correlations in the wave function of the ground state. 
Therefore, we anticipate that the value of A g we found is robust: it is not an artifact of the 
Maki-Zotos trial function, but it should be the same (with logarithmic accuracy) for any 
state homogeneous on the macroscopic scale r > R c . 



VI. SPIN EXCITATIONS 

In this section we analyze the simplest excitations of the PFLL at Av = 1. This problem 
is similar to the one studied in Ref. |14|. The only difference is that the interaction potential 
is renormalized due to the polarizability of the other Landau levels. Our goal is to study 
how this renormalization affects the spin excitations. 

The energy of the fully spin polarized ground statfl \M = Md?, S z = M*/2) at Av = 1 

is 



n - m$ U* + 



A ea; + gohuj c 



(6.1) 



where A ex is defined by Eq. ( |4.4| ). The eigenstates of the Hamiltonian with an extra electron 
or hole 



M 4 



\h) = il>N,k,t 



Md>, 



Mi 



have the energies 



E l = E p ~ I 1 * + \9&u c , 
E } = E p + fj,* + A ex + y o huj c 



(6.2) 
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respectively. Equations ( |6.2|) enable us to relate the width of the spin gap A s = E± +E^ — 2E P 
to the value of A ex , i.e. A s = A ex + g hu> c . Substituting A s in Eq. fl2.9|) and using Eq. Q4.4|) , 
we obtain the result ( p . 1 0| ) for the effective ^-factor. 

It is worth noticing that the effective ^-factor is independent on N in the limit of a weak 
magnetic field due to the Thomas-Fermi screening^. (Without screening, this factor would 
logarithmically diverge with N —>■ oo.) 

The exchange enhancement is maximal if the filling factor is odd. When the both spin 
sublevels are either empty or completely filled iy is even), the spin splitting is determined 



emp 



by the bare g - factor 

Let us turn to the consideration of neutral excitations - spin wavesBS which at Au = 1 
are described by the wave functions: 



|sw,k> = 5> 



(6.3) 



Wave functions ( |6.3D are eigenstates of the Hamiltonian (3.27) with the energies E p + Esw(k) 
where the energy of the spin wave Esw{k) is given by 



Eswik) = g c shu c - 



(6.4) 



d 2 q 2ne 2 1 
(2tt) 2 nq 7{qj 



2l)2' 



q 2 t 



-i 2 



-<j 2 ^ 2 /2+ikq^ 2 



For small wave vectors, kR c 1, the spin wave has a quadratic dispersion relation which is 
typical for ferromagnets: 



e 2 R. 



Esw(k) = g huj c + '—^-k 2 



71 K 



(6.5) 



At k — > the energy of the spin wave is determined by the bare ^-factor in agreement with 
Larmor's theorem. It is worth mentioning that the long- wavelength domain of the spin wave 
spectrum is controlled by the bare Coulomb interaction. This is because the spatial scale 
k£ 2 important for the spectrum, see Eq. ( |6.4D , is much smaller than the screening radius a#. 
In the region R~ l <C k -C /cf the energy of the spin wave reveals oscillatory behavior 



Esw(k) = g cS hu} c + 



(ke\ _ sin(2kR c 



a B 



2kR r 



(6.6) 



Here C(x) is a smooth function with asymptotes 



C(x) 



x < 1 



Spin waves with extremely large wave vectors, k 3> kp, correspond to almost independent 
electron and hole, see Eq. (|6.2|). In this limit, the energy Esw approaches A s : 



A s - E sw (k) = U(k£ 2 ) 



2 \ < 



flUlc 

2N 



ln Nr* 



(6.7) 
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VII. CONCLUSION 



In this paper we constructed a theory describing low-energy excitations in a 2D electron 
liquid in a weak magnetic field. We have shown that all the excitations with energy smaller 
than cyclotron energy huj c can be described by the effective Hamiltonian acting in the Fock 
space of the partially filled Landau level only. Starting from first principles, we obtained 
the explicit form of this Hamiltonian by integrating out all the other degrees of freedom. 

Armed with this effective theory, we have been able to make important predictions: 

1) We found that the thermodynamic density of states is negative for non- integer filling 
factors. For a broad range of filling factors, the value of the thermodynamic density of 
states was found to be independent on the value of magnetic field, see Eq. ( |4.6|) . This effect 
may be revealed in magnetocapacitance measurements 2 ^ in the weak magnetic field regime. 

2) The tunnelling density of states was shown to have a gap at the Fermi level. For a broad 
range of filling factors, the width of the gap, see Eq. ( J2.8p , was shown to be consistent with 
that predicted by hydrodynamic approach of Ref. []. The gap can be observed in tunnelling 
experiments!! in a weak magnetic field. Evidence of suppression of the tunnelling density of 
states for the filling factors v ^5 9 was reported recently by Turner et. a/.EI. This suppression 
may be associated with the gap predicted by our theory. However, the observed width of 
the gap is twice larger than that given by Eq. ( j2.8| ). 

3) The exchange enhancement of the effective g-factor remains strong in the weak field 
regime, and at all odd filling factors it takes a universal value, see Eq. fl2.10| ). We found 
that the energy scale of charge excitations is parametrically smaller than the energy scale 
for spin excitations. This is qualitatively different from the situation at low filling factors, 
v 1, where both these excitations are characterized by the same energy scale e 2 fnl. 

None of the aforementioned effects could be obtained in the framework of the Landau 
quantization of the spectra of quasiparticles in the conventional Fermi liquid theory. This 
theory describes adequately the excitations with energies larger than fuo c . Thus, our de- 
scription of the low-energy properties of interacting 2D electron gas is complementary to 
the Fermi liquid picture. 
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APPENDIX A: EVALUATION OF 

Integration in Eq. ( |3.21| ) immediately yields 



n(g,*) = -^£ £ 

n m=0 n 2 =N L 



(-l)(*w-»i)(n2 -7li) 
{uj/uJc) 2 + (n 2 - nx) 2 
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j 712— n 
L n 1 



2 02- 



jn 1 -n 2 ( 1 ^ 

K 2 I — 



-q 2 P/2 



(Al) 



where L^(x) is the Laguerre polynomial^!. 

Further simplification is possible for q <C kp. Under this condition, we can use the 
asymptotic expression for the Laguerre polynomial^ 



L™(x) 



in 



+ m)\e x / 2 J m ^2x(2n + m + 1 



n\[x(n + m/2 + l/2)] 



m/2 



(A2) 



where ^T m (x) is the Bessel function. This expression is applicable if x <C n. Using Eq. (|A2| 
and introducing a new index j = ri2 — n\ we approximate expression ( |AT] ) as: 



U(q,uj) = - 

N-l 
ni=max(iV— j',0) 



2m ~ 



J, 



(uj/uj c ) 2 +f 



2n x + j + 1 



(A3) 



The terms giving the main contribution to the sum in Eq. ( |A3| ) are characterized by value 
j ~ q£\/~N <C iV. This enables us to approximate 2n + j + 1 « 2iV and perform the 
summation over n\. This yields^ 



, . 2m ~ 

n(g,o;) = E 



7T ^ (w/^ c ) 2 + J 



2 L^J 



(A4) 



Equation (|A3|) can be transformed into an integral form: substituting the identity 



into Eq. 



[X(^)] 2 = / — cos ny Jo ( 2x sin ~ 
Jo i V 2/ 

) and performing the summation over j, one easily obtains: 

y\ ^COsh^u" 



m 



71 



dyj 2qR c cos 



sinh7r— 



(A5) 



(A6) 



Equation (|A6|) enables one to obtain the asymptotic form of the polarization operator in 
different regimes. In the static limit, to <C u c , one finds from Eqs. (|A6|) and (|A5|) : 



U(q,uj) = - 



m 



71 



1 - J 2 (gi? c ) + O 



10 



tot. 



(A7a) 



For the high-frequency domain, ^> u c , only a small vicinity of the point y = it contributes 
to the integral in Eq. (|A6|), and H(q, uS) coincides with the result for 2D electron gas in the 
absence of the magnetic fieldS 



IL(q,u) = 

7T 



\UJ\ 



u 2 + q 2 vp 



(A7b) 
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In the hydrodynamic limit, qR c <C 1, one can expand the Bessel function in Eq. (|A6|) in a 
Taylor series, which yields 



U{q,u) = - — 

Z7T 



q 2 v 2 F 
UJ 2 + uj\ 



+ 0{q 4 Rt) 



(A7c) 



In the opposite limit, R c 1 <C g <C &f, the main contribution to the integral in Eq. ( |A6 
comes from the vicinities of points y = and y = ix. The calculation gives 



Yl(q,u) w - 



??? 



7T 



COth f 7T — 



^Juj 2 ~+~q 2 ~v§ 



LO 



sin 2gi? c 



S( 9 ,a;) = l + 



gt>i? ) sinh 7r^- ' 

q 2 VpL0 2 (4a; 2 — g 2 f£0 



,2\3 



■ (to 2 + q 2 v 2 F ) 



(A7d) 



APPENDIX B: PARAMETERS OF THE EFFECTIVE HAMILTONIAN 



In this Appendix we present the details of the derivation of the effective Hamiltonian 

We start from Eq. ( p.22j ). We notice that operators F\ and F^, see Eqs. ( |3.11c|) and 
( p.lld|) , are integrals of the Matsubara operators over the time interval [0, j3 }. Let us divide 
this time interval into M small intervals [r m _i,r m ] , m 



1,2. . .M - 1, where t„ 



mAr 



and At = (3q/M. We imply that At must be smaller than the characteristic time of the low- 
energy dynamics, but much larger than uj' 1 . Then, the chronologically ordered exponent in 
Eq. (|3.22|) can be factorized as 



T T {e-*-**}= f[T T \e^ m) -^ m) \ 

m=l ^ > 



(Bl) 



where operators F^ and F^ have the form 



(m) 



,(m) 



-1 Ft. 



:b2) 



Definitions (|B2|) differ from the corresponding definitions ( 3.11c ) and (|3.11d| ) by the con- 
straint on the time domain of integration. We do not need to impose an additional con- 
straint on the time component of £2 in the second of Eqs. (|B2|) because Go decays rapidly 
at |£? - £2! > 1M, see Eq. fl37TJ). 

Now, we substitute Eq. (|BT| ) into Eq. ( |3.22|) and perform the integration over the field 
cf). It yields 
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A(/3b) = e -*(no+AQ + ^) -q / Tr L-f^-f^XX , (B3) 

m=l ^ * ' 

where average (J 7 ) of an arbitrary functional T {(j)} over the field cf) is defined as 

(JF) = e PoAn Jv<f) ^{<j)}x 

exp Qy rf 3 ^ | rf 3 60(^i)0(6)^(6 -6)) • (B4) 
The normalization factor e /3 ° An describes the fluctuations of field <p around the saddle point: 
e~ PoAn = jv^e^Jd^Jd^m^S^-i')) ■ (B5) 



An explicit expression for Ail is presented in Subsection |B 1| . 

Because At is small, we can expand every factor in Eq. (|B3| ) up to the first order in At: 

T T je-^-^j^ « l _ AtHj {* N (r n ), ^(r n )} . (B6) 

This relation defines the operator Hi. Let us notice that this operator actually does not 
depend on r n because it contains the same number of creation and annihilation operators 
taken at the same moment of time. Having this in mind, we substitute Eq. ( p6|) into Eq. ( |B3|) 
and obtain 

A(/3b) = e -^ +Af ^) _ ^ Hl ) . (B7) 

Finally, we take the limit M — > oo in Eq. ( |B"7| ) and substitute the result in Eq. ( |3.6| ). It 
yields 

H eff = Q + AQ + H£ + H T . (B8) 

Now, one has to perform the actual calculation of the average in the left hand side of 
Eq. (|B6| ) in order to find the operator Hj. It can be done by using perturbation theory. The 
leading terms in the small parameters r s and 1/N are 

T r {e-^ ) -^ ) })^l-(F 1 W + F 2 (i 

+ 1 -(T T (Ft ) + Fi n) ) 2 ). (B9) 

To find the operator Hj, only terms linear in At should be retained. 

Calculation of the right-hand side of Eq. (|B~9| ) is carried out with the help of Eqs. ( |B~2"D 
and of the correlation functions of the field <p: 

(0) = 0, (BIO) 

(0(r,r)0(O,O)) = -^(r)- 
K\r\ 



(2tt) 3 ^ nq ) e(q,u) 
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The latter expression follows from Eq. ([B4] ) and ( |3.19|) . In the second of Eqs. ( P10|) , we 
explicitly separated a term which is singular at r = 0. 

The averages that give the leading contribution to the effective Hamiltonian are 



■<n) 



1 , Tt { p [n) 



At (/i ex . + /x c ) / d 2 rp N (r) 1 

= -Ar (hHI + ^ p J d 2 rp N (r] 



(Bll) 



where the operator i^ 6 // is defined by Eq. (|3.23|) . 

In Eqs. ( |B11[ ), \x ex coincides with the well-known exchange correction to the chemical 
potential: 



2 

fie* = -2irf ( d 2 r—G (r, 0; 0, e -> +0) P N (0, r) 



(B12) 



where Go is defined by Eq. ( |3.16| ), and /i c , p: p are the correlation corrections to the chemical 
potential found in the random phase approximation: 



dud 2 q f2ne 2 \ U(q,uj) 



(2tt) 3 V 

x (_ 1 )(AT-n) jL Ar-n| 



'q 2 £ 2 ' 



E 



u c {N - n) 



■n-N / Q 



2 1)2 



-q 2 £ 2 /2 



fi p = Tr£ 2 J d 2 rP N (0, r) P N (r, 0) (tf (r) - ^ . 



(B13) 
(B14) 



When deriving Eqs. ( gig ), ( |BT3| ), we use Eqs. (^TJ and (|3TT6| ). 

Finally, substitution of Eq. ( pil| ) into Eq. ( P9| ) enables us to find the operator Hi. 
Comparing Eq. (|B8| ) and Eq. (|3.27|) , we obtain the expressions for the parameters of the 
Hamiltonian H e ff. 



n = n + An, 

p, - Nhu c - fi ex - fi c - fip. 



(B15) 
(B16) 



The following subsections are devoted to the explicit calculation of the parameters of the 
effective Hamiltonian in terms of the filling factor and interaction strength. We will assume 
that condition ( 4.18|) is met. 



1. Calculation of AQ. 



Equations (Ml 



D and (|3~T9l ) yield 

AQ 1 r dud 2 q 
1J ~JZkY 



L X Ly 



In (s(q,uj)) . 



(B17) 



With the help of Eqs. ( |3.20|) , ( |3.21| ) and ( |3.14|) , one can transform Eq. (pi7|) into a more 
convenient form: 
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AQ = AQ ex + AQ 

An 



L X Ly 

Att c 

L X Ly 



d 2 r—G (r } 0;0,e) G (0, 0; r, e) 



1 /"i /-dwt^g /2vre 2 \ 2 



"27o ^ (2vr^ 



all 2 (g, cu) 



i-a¥n(p) 



(B18) 
(B19) 

(B20) 



Expression ([B19| ) corresponds to the first order exchange correction to the ground state 
energy, and e — > +0. Equation 



is the correlation correction to the ground state 
energy equivalent to the sum of ring diagrams^. 

First, we evaluate the exchange energy ( |B19| ). Using the explicit expression for the 



Green's function ( |3.16| ), the property of the Landau level wave functions, 



Pn(r 1 ,T 2 ) = X)^n,fe( r 2)^n,fe(ri) = 
k 

1 D i(yi-y2)(^+ x 2)/2£ 2 -\r 2 -r 1 \ 2 /4e 2 £ 



2ir£ 2 



|£a - ri| 

2£ 2 



and the identity for the Laguerre polynomials, 



(B21) 



i=0 



L„ (x) 



(B22) 



we obtain 



ATLex 

L x L y 



2 

dxe~ x 12 



(2tt£ 2 ) 7 

Integration in Eq. ( |B23|) can be easily performed, which yields for N 3> 1: 



(B23) 



L X Ly 



KTT 2 £ 3 



(2N) 



3/2 



lnA^ 



16(2iV)V2 \W/ 2 



(B24) 



The first term in brackets coincides with the exchange energy of 2D electron gas at zero 
magnetic field, and the second term appears due to the confinement of the electron wave 
functions by the magnetic field. The logarithmical factor IniV ~ In (kpR c ) arises due to 
the integration over relatively large spatial scales hp 1 <C r ^$ R c in Eq. ( |B19| ). On these 



scales, the screened interaction potential is significantly smaller than the bare potential 
appearing in Eq. ( |B19| ). Therefore, it is plausible to anticipate that the correlation term 
of the thermodynamic potential, which accounts for the screening effect, should partially 
compensate the large logarithmic factor in Eq. (|B24j ) , and lead to the replacement R c — >• Ob 
in the argument of the logarithm. The calculation of the correlation energy ( |B20| ) explicitly 
demonstrates this. 

In zero magnetic field the contribution of the correlation energy^ is only of the order 
of 0(r 2 ). The situation changes in the magnetic field, where a contribution proportional 
to r s ln(Nr s ) appears. This term arises mainly from integration over the domain of wave 
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vectors i? c 1 C g <C a B l in Eq. (|B20|) , where we can use the asymptotic expression ( |A7d|) 
for n(g, uj). It yields with the logarithmic accuracy 



Aa 



L x L y K lV2N (16tt£ 2 ) 



\n(Nr s ) + 0(r 2 s ). 



(B25) 



Finally, with the help of Eqs. (|B24j ) and (|B25| ), we obtain for the thermodynamic poten- 
tial (|BT5| ), 



_ LxLy_ 



' iV(jV-l) 



ttk£ 



16(2JV) 



1/2 



In (r; 1 ) 



(B26) 



This result contradicts Ref. |32| where the contributions (|B19|) and ( [B20|) were evaluated and 
a correction oc r^ 4 was obtained instead of rJ InrJ. 



2. Calculation of //. 

We start from the calculation of the exchange correction to the chemical potential. With 
the help of Eqs. (|37T6| ), (|B21| ) and (|B22|) we transform Eq. (|BT^ ) to the form 



^ /" cte-^Ljr-i ( ^ ) M 



a; 



(B27) 



After integration in Eq. ( |B27| ) we obtain for N ^> 1 with logarithmic accuracy: 

2 „. 1 In 1 ^ 



^ ~ Ki U (2iV)1/2 2tt(2A0 1 /2 



In AT 



(B28) 



The correlation shift of the chemical potential, Eq. ( |B14[ ), is calculated with the help of 
Eqs. ( |B21|) and ( |3.25|) . With the logarithmic accuracy, the calculation gives 



k£ 27r(2iV) 1 / 



^ In (iVr s ). 



(B29) 



It can be checked by an explicit calculation that the other correlation correction to the 
chemical potential, /i c , see Eq. (|B13| ), contains an additional small factor in comparison with 
Eq. ( |B29| ): fi c ~ fJ> p /N, and, therefore, it can be neglected. Finally, with the logarithmic 
accuracy, we have for \f appearing in Eq. (|Bl!) the result: 



e l (2 ,_ 1/2 



/i - Nhuo c + — - (2JV) 

KZ \ 71 



In (r„ 



2vr(2iV) 1 /2 y • 



(B30) 
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FIGURES 



FIG. 1. The renormalized pair interaction U(r) given by Eq. ( 3.25| ) as the function of the 
inter-electron distance (solid line). Dashed line is the bare Coulomb potential. The asymptotic 
behavior of the renormalized potential is given by Eqs. (|2.5[) and (2.6) 



FIG. 2. The thermodynamic density of states as a function of the filling factor of the partially 
filled Landau level in the domain < v < 1. The asymptotic behavior of the thermodynamic 
density of states is given by Eqs. (Op, (KW), and (|42lD 



FIG. 3. (a): The spectral density as a function of energy, (b): The dependence of the 
characteristic energy scales on filling factor of the partially filled Landau level Av. The asymptotic 
behavior of the width of the pseudogap, A g is given by Eq. ( |5.12j ). The width of the peak in the 
one-electron spectral density A p is determined by the variation of the interstitial energy and it is 
smaller than Twu c /2N at all filling factors. 
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